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Abstract 



We endow a recently devised algorithm for generating exact eigensolutions of 
large matrices with an importance sampling, which is in control of the extent 
and accuracy of the truncation of their dimensions. We made several tests 
on typical nuclei using a correlated basis obtained from partitioning the shell 
model space. The sampling so implemented allows not only for a substantial 
reduction of the shell model space but also for an extrapolation to exact 
eigenvalues and E2 strengths. 
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I. INTRODUCTION 



New effective methods for solving exactly the shell model (SM) eigenvalue problem in 
complex nuclei have been developed in recent years. A notable one is the Monte Carlo 
technique [1] for generating exact SM ground states [2-4]. Its application, however, has 
been rather restricted because of the well known minus-sign problem. 

Alternative methods face directly the diagonalization of the SM Hamiltonian by resorting 
to algorithms like Lanczos [5] and Davidson [6] . The critical points of direct diagonalization 
methods are the amount of memory needed and the time spent in the diagonalization process. 

In order to overcome these limitations there have been attempts to combine the stochastic 
methodology with the standard diagonalization approach [7-12]. In [7] the use of Gaussian- 
like single-particle basis states with random oscillator frequencies was suggested. The au- 
thors of [8-10] proposed a stochastic truncation of the SM matrices and suggested for the 
energies so obtained an exponential extrapolation law to the exact eigenvalues. In the so 
call quantum Monte Carlo diagonalization method (QMCD) [11,12] the reduction of the 
dimension of the SM space is achieved by using the auxiliary field Monte Carlo technique to 
select the relevant basis states. Though quite successful, the latter method is not free from 
problems. Indeed, one has to deal with the redundancy of the basis states, which may slow 
considerably the convergence of the procedure, as well as with the problem of restoring the 
symmetries generally broken in stochastic approaches. 

In a recent paper [13], we have developed an iterative algorithm for determining a selected 
set of eigenvectors of a large matrix which is faster than the other, currently adopted, 
algorithms, including Lanczos, and extremely simple to be implemented. It is, moreover, 
robust, yielding always ghost-free stable solutions. 

Like the other methods, however, it requires the storage of at least one eigenvector, 
which exceeds the limits of modern computers in many complex systems. The present 
paper deals with the problem of overcoming this limitation. To this purpose, we endow the 
algorithm with an importance sampling for reducing the sizes of the matrix. The sampling 
is closely linked to the iterative algorithm and is in full control of the accuracy of the 
eigensolutions. This will emerge clearly from the results of the exhaustive tests presented 
and discussed in this paper. In fact, we adopted the importance sampling to solve the shell 
model problem for three typical nuclei, the semi-magic 108 Sn, the N = Z doubly even 48 Cr, 
and the N ^ Z odd 133 Xe. In order to enhance its efficiency we used a correlated basis 
obtained from partitioning the shell model space according to a method developed in Ref. 
[14] . We will show that the method not only allows for a drastic truncation of the SM space 
but yields naturally extrapolation laws to exact eigenvalues and eigenvectors as well as to 
the occupation numbers and the E2 transition probabilities. 



II. THE ALGORITHM 

Let us first give a brief outline of the algorithm [13]. For the sake of simplicity, we 
consider here a symmetric matrix 

A = {a ij } = {{i\A\j)} (2.1) 
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representing a self-adjoint operator A in an orthonormal basis {| 1), | 2), . . . , | N)}. The 
algorithm goes through several iteration loops. The first loop consists of the following steps: 
la) Diagonalize the two-dimensional matrix (a^) (i,j=l,2), 
lb) select the lowest eigenvalue A 2 and the corresponding eigenvector 

I0 2 )=cS 2) |1) + 4 2) |2), (2.2) 

lc) for j — 3, . . . , N, diagonalize the matrix 

where bj = (4>j-i I A \ j) and select the lowest eigenvalue Xj together with the corresponding 
eigenvector | <f>j). This zero approximation loop yields the approximate eigenvalue and 
eigenvector 

N 

E^ = \ N , |^ (1) )=l0 J v)=E4 7V) M)- (2-4) 

i=i 

With these new entries we start an iterative procedure which goes through n = 2, 3, . . . 
refinement loops, consisting of the same steps with the following modification. At each step 
j — 1, 2, . . . , N of the n-th loop (n > 1) we have to solve an eigenvalue problem of general 
form, since the states | <f>j-i) and | j) are no longer orthogonal. The eigenvalue = A^ 
and eigenvector | ip^) =| obtained after the n-th loop are proven to converge to the 
exact eigenvalue E and eigenvector | ip) respectively [13]. 

The algorithm has been shown to be completely equivalent to the method of optimal 
relaxation [15] and has therefore a variational foundation. 

Because of its matrix formulation, however, it can be generalized with minimal changes 
so as to generate at once an arbitrary number v of eigensolutions. Indeed, the first loop goes 
through the following steps: 

la) Start with m{> v) basis vectors and diagonalize the m-dimensional principal submatrix 
Rj} = l,m), 

lb) select the v lowest eigenvalues Ai, A 2 , . . . , A^ and the corresponding eigenvectors | | 

02), • • • , | 4>v)- 

For k — 2, 3, . . . , fcjv, where kN are the steps necessary to exhaust the whole iV-dimensional 
matrix, 

lc) diagonalize the matrix 

where is a f-dimensional diagonal matrix whose non-zero entries are the eigenvalues 
\[ k 1 \ A 2 fe x \ . . . , A^ fc_1 ), Ak = {aij} = (k — l)p+l, . . . , kp) is a p-dimensional submatrix, 
B k and its transpose are matrices composed of the matrix elements b\® = (0^ ^ | A \ j) (i 
= 1, j = (k-l)p+l,...,kp). 

Id) Select the v lowest eigenvalues and the corresponding eigenvectors 

| # } >,| #>,...,! 0( fc )>. 
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Once the basis is exhausted, the process yields v approximate eigenvalues and eigen- 
vectors, which are the new entries for a new iteration. This goes through the same steps 
with one essential modification. Each loop, in fact, can be viewed as the solution of the 
eigenproblem for the restriction A\ s of the operator A to a subspace defined by the span of 

the set of vectors | i p ) = {(f)[ k \ . . . , ($\ \ (k — l)p + 1), . . . , | kp)}. Since the vectors <f>f^ 
are linear superpositions of all the basis vectors | j) we started with, the basis obtained 
after the first approximation loop is no longer orthonormal, just as in the one-dimensional 
eigenspace, and may be even redundant. We have therefore to solve an eigenvalue problem 
of general form. This is done most effectively through a pivotal Choleski decomposition of 
the overlap matrix (i p \ i' ) [16]. As outlined in a forthcoming paper, this procedure is very 
easy to be implemented and extremely fast, specially in view of the very limited dimensions 
of the matrices coming into play. With this modifications, the subsequent iteration loops 
proceed as the first one. In this way we generate a sequence of v vectors ^S fc) ,...,^ fc ). The 
restriction of the operator A to these sets defines a sequence of diagonal matrices, whose 
non zero elements are the current eigenvalues X^\...,X^\ with decreasing norms. This 
monotonic sequence is certainly bounded from below and therefore convergent. 



III. IMPORTANCE SAMPLING 

The just outlined algorithm, though of simple implementation, requires the storage of at 
least one eigenvector. Since for many complex systems the dimensions of the Hamiltonian 
matrix become prohibitively large, one must rely on some importance sampling which allows 
for a truncation of the space by selecting only the basis states relevant to the exact eigenso- 
lutions. A notable example is the stochastic diagonalization method [6], which samples the 
basis states relevant to the ground state through a combination of plane (Jacobi) rotations 
and matrix inflation. 

A similar, but more efficient, strategy can be implemented in the framework of our di- 
agonalization process. Exploiting the fact that the algorithm yields quite accurate solutions 
already in the first approximation loop, we can devise a sampling which makes use of the 
first loop only. This is to be accordingly modified and goes through the following steps: 

la) Turn the v- dimensional principal submatrix {a^} = l,v) into the diagonal form 
A„with eigenvalues Ai, A2, . . . , A„. 

lb) For j = v + 1, . . . , N, diagonalize the v + 1-dimensional matrix 

A = ( £ S > ) , (3.1) 

where bj = {b lj} b 2j , •, b vj }. 

lc) Select the lowest v eigenvalues A^, (i = l,v) and accept the new state only if 

El^-A,|>6 (3.2) 

i=l,v 

Otherwise restart from point lb) with a new j. 

We can avoid such a time consuming sampling procedure by resorting to an alternative, 
though completely equivalent, route based on the method developed for deriving an exact 
non perturbative shell model Hamiltonian [17]. We carry the similarity transformation 
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A' = Q^AQ (3.3) 

where 

Here, I v is the f-dimensional unit matrix and uo a f-dimensional vector. 
The transformed matrix has the following structure 

\ h j "l, - • h J J 

where 

b'j = — (uj ■ bj)cu — ljA v + cijjLd + bj. (3.6) 
We now impose the decoupling condition 

b'j = — (uj ■ bj)uj — ljA v + cijjLJ + bj = 0. (3.7) 

Once we find the solution u, the matrix element ajj — uo -bj becomes an eigenvalue of the 
matrix A' and, therefore, of matrix A. Right multiplying Eq. (3.7) by bj, we obtain the 
dispersion relation 

i=i, v a jj — \i — to ■ bj 

This admits v + 1 solutions, corresponding to the v + 1 eigenvalues of A. In correspondence 
of the lowest solution (Co ■ bj) min , we get the maximum eigenvalue 

X 'max = a ii - • bj)min- (3.9) 

The eigenvalues = 1, . . . ,v separate at least in weak sense the new eigenvalues X'^i — 
1, . . . , v + 1 [18], namely 



Since 



A; < Ax < A' 2 < A 2 < . . . < A: < A, < X' max . (3.10) 



£ K = £ A- + ajj - (w • fejOmin = ajj + X i ( 3 - n ) 

i=l,v+l i=l,v i=l,v 



we have 



E (A- - Ai) = (u ■ bj) min . (3.12) 

1=1,1; 

Condition (3.2) is therefore equivalent to 
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(to ■ bj) min > e. 



(3.13) 



The just outlined sampling procedure requires only the solution of the dispersion equation 
(3.8), which is of the type 

f(z) = z (3.14) 

and fulfils the condition 

l-df(z)/dz>0. (3.15) 

We can therefore easily solve it by using the Newton method of derivative. This alternative 
sampling procedure is not only rigorous but also more economical. It avoids, in fact, the 
(N — v)— fold iterated diagonalization of v + 1 dimensional matrices. 



IV. MULTIPARTITIONING METHOD 

The extent of truncation induced by the sampling is maximal when the eigenvectors are 
highly localized. This condition is fulfilled in most physical problems. Even when this is 
not the case, we can approach the above condition by using a correlated basis obtained by 
a multipartitioning method [14]. This goes through the following prescriptions: 
i) Partition the shell model space for N valence nucleons into orthogonal subspaces, Pi and 
P 2 according to 

P = Pi + P 2 , (4.1) 



ii) distribute Ni and N 2 nucleons (Ni + N 2 

ways, 

iii) decompose the Hamiltonian H into 



= N) among these subspaces in all possible 



H = H 1 + H 2 + H 12 , (4.2) 

iv) solve the eigenvalue equations 

Hi | a^) = E ai | ai Ni) (4.3) 

obtaining the eigenstates | a±Ni) and | a 2 N 2 ) of Hi and H 2 respectively in Pi and P 2 . Once 
this is done, it is possible to replace the standard shell model basis with one composed of 
the states 

| a N) =| aiA^ia 2 A^ 2 ). (4.4) 

We use the above basis to diagonalize the residual term Hi 2 of the SM Hamiltonian. The 
new basis is in general highly correlated and, therefore, highly localized in the Fock space, 
a feature which enhances considerably the efficiency of the method. 



6 



V. SELECTIVE NUMERICAL TESTS 



We applied the sampling algorithm to the semi-magic 108 Sn, the N=Z even-even 48 Cr 
and the N > Z odd-even 133 Xe. The model spaces are: 

1) P = {2d5/2, lgl/2, 2d3/2, 3sl/2, l/ill/2} for the 8 valence neutrons of 108 Sn and for the 
4 valence protons and 3 valence neutron holes of 133 Xe, 

2) P = {1/7/2, 1/5/2, 2p3/2, 2pl/2} for the 4 valence protons and neutrons of 48 Cr. 

We adopted a realistic effective interaction deduced from the Bonn-A potential [19] for 
108 Sn and 133 Xe, and used the KB3 interaction [20] for 48 Cr. For 108 Sn we used the single 
particle (s.p.) energies (in MeV) e 2d5/2 = 0, e lg7/2 = 0.2, e 3s i /2 = 2.2, e 2d3/2 = 2.3, e lhll/2 = 
2.9. Apart from the first two currently in use, the other energies were deduced from a fit 
on m Sn for {2d3/2, 3sl/2} and from the level 9~ in 106 Sn for l/ill/2. We used the same 
set of energies for protons in 133 Xe, while for neutron-holes we used the slightly different 
values e 2d5/2 = 0, e lg7 / 2 = 0.2, e 3sl/2 = 1.72, e 2d3/2 = 1.88, e mi/2 = 2.7. The different 
set accounts effectively, at least in part, for the significant asymmetry between proton and 
neutron numbers which yield different contributions from the three- and more-body forces, 
not included in the effective interaction. As for 48 Cr, we used the same energies adopted in 
Ref. [20], namely e lf7/2 = 0, e 2p3/2 = 2.0, e 2p i/ 2 = 4.0, e lf5/2 = 6.5. 

We partitioned the shell model space for 108 Sn according the following prescriptions: 

/ P\ = {2d5/2, lgl/2} 
P = {2^5/2,1^7/2, 2d3/2, 3sl/2,l/ill/2} + (5.1) 

\P 2 = {2d3/2,3sl/2,l/ill/2}. 

This partition is dictated by the large energy gap (~ 2 MeV) between the two corresponding 
sets of single particle energies. 

For 48 Cr and 133 Xe, we simply decompose the space into a proton and neutron subspace 

P = P P + P n . (5.2) 

We adopted then the multipartitioning method [14] to generate a new correlated basis 

\j)=\aN) (5.3) 

and used this new basis to implement the importance sampling. 



A. Eigenvalues 



For illustrative purposes we discuss only few of the lowest states of the nuclei under 
investigation. As shown in the plot of Fig. 1, the sampling parameter e varies very smoothly 
with the dimensions n of the reduced matrices. In these, as in all other states considered, it 
scales according to 



e = b— exp 

n 2 



N' 
— c — 
n . 



(5.4) 



In virtue of this law, the increment of the dimensions of the matrix is modest for large 
values of e, but grows dramatically as e gets smaller and smaller. This behavior reflects 
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the density of the unperturbed levels which is very low at low energy and raises steeply 
around a centroid at high energy. This is shown in the upper panel of Fig. 2. It follows 
that the running sum of the unperturbed basis states grows very slowly at low energy, 
then raises steeply toward its saturating full N value in a relatively small energy interval 
(middle panel of Fig. 2). It is important to point out that this is the distribution of our 
unperturbed correlated states defined by Eq. (4.4). The partitioning of the shell model 
space is responsible for the squeezing of their energies around a centroid. Had we adopted 
the standard SM basis, the state distribution would have resulted more spread out and the 
running sum would have grown discontinuously through several steps in correspondence of 
each sub-shell (lower panel of Fig. 2). 

Figs. 3-5 plot the eigenvalues versus the dimensions n of the matrices resulting from 
decreasing values of the sampling parameter e for some low- lying states of 108 Sn, 48 Cr, and 
133 Xe, respectively. 

In all nuclei and for all states, the eigenvalues decrease monotonically and smoothly with 
n. The only meaningful exceptions are represented by the curves of the first excited J 71 " = + 
and J n = 2 + and few other, similarly behaving, high-lying states of 48 Cr. These undergo 
a jump from an upper to a lower curve at some small value of the sampling parameter 
e, a signal of energy crossing. Even in these cases, however, the subsequent behavior of 
the energies is smooth as for the other states. It follows that, in all cases, starting from a 
sufficiently small e, the energies scale with the dimensions n according to the law 



N 

E = E + b — exp 

n 



N 
-c — 

n 



(5.5) 



where b, c, and E are constants specific of each state and the full dimension N provides the 
scale. This fit allows for an extrapolation to asymptotic eigenvalues which differ from the 
exact ones in the second or third decimal digit. This is explicitly proven for 108 Sn and 48 Cr 
(Table I). The reliability of the extrapolation in the case of 133 Xe is inferred from the rapid 
convergence of the iterative procedure, as clearly illustrated in Fig. 5. Indeed, the curves 
reach a plateau of practically constant energies starting from an n value which is smaller 
than the full dimension N by more than one order of magnitude in 108 Sn, 48 Cr, and by more 
than two in 133 Xe. 

Our exponential extrapolation law is somewhat different from the one proposed in Refs. 
[9,10]. On the other hand, it is valid for all states and nuclei examined and follows directly 
from the sampling, as it can be inferred from the following heuristic argument. 

Let us consider the simplest case of one-dimensional eigenspace (v — 1). From Eqs. (3.8) 
and (3.12), we can write the corresponding sampling prescription as 

AA = E AA, = E(A; - Aj) = - E j~ AA~ (5 ' 6) 

j j j " j j ~ A i - AA i 

Expanding A A,, we get a solution whose leading term is 

AA, = - (5.7) 

/,From the plots it is clear that the extrapolation law holds for an energy range of 1-2 MeV in 
proximity of the exact eigenvalue. It accounts therefore for small contributions coming from 
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a small fraction of the basis states in the J2j m the high density region around the peak, 
as shown in Fig. 2. Since, in this range, the energies Ej of our correlated basis states are 
closely packed and approach the E n value, we can put (ajj — A) oc n in first approximation. 
It remains to analyze the matrix elements b\y These are given by 

(0, | v | j) 2 (5.8) 

where 

|&> = E<i|i). (5-9) 
i=i 

For the lowest eigenvalue, the dominant q components of | 4>j) are the low-energy ones 
with small % values. It follows that the probability, b\^ that the interaction couples | <f>j) to 
| j) is small and random for j < n and vanishes for j > n, as prescribed by the sampling 
criterion. This implies that the dimension n represents the range of the allowed events. We 
can therefore put b\j oc exp(—j/n), where j is a label representative of the small fraction of 
j terms of the sum J2j- We used the factor N to fix the scale. The scaling law (5.4) for e 
follows from the one for the energy E since e is essentially the derivative of E with respect 
to n. 

We have also compared the sampling with the energy truncation of the Hamiltonian 
matrix. As shown in Fig. 6, the sampling is obviously more effective and accurate. The 
two procedures, however, tend to become equivalent for a value of n which, though large, 
is still much smaller than the full dimension N. It is to be pointed out, however, that, in 
our case, the effectiveness of the energy truncation is due to a great extent to the use of the 
correlated basis obtained through the multi-partitioning method. 



B. Eigenvectors and E2 transitions 

An accuracy of the same quality is reached for the eigenstates of the n-dimensional 
truncated Hamiltonian matrix 

n 
i=l 

where | i) are the correlated basis states obtained by the partitioning method. 

In Fig. 7 we give the the overlap of ip n with the exact eigenvector ip for the first five 
J 71 " = 2 + states of 108 Sn and J 71 " = + of 48 Cr. The convergence to unity is fast for all five 
states, even if, for some of them, the overlap is very small at small n. Small fluctuations 
are noticeable at small n values. They reflect the interference between the components 
of different wave functions in correspondence of partial energy crossings. The above two 
features represent a further proof of the robustness of the iterative algorithm. 

To complete the analysis we study the behavior versus 1/e of the strengths of the El 
transitions between some low-lying states in 108 Sn and 48 Cr ( Fig. 8 ), as well as in 133 Xe 
(Fig. 9). In all cases, the strengths reach soon a plateau and, then, undergo very small 
variations, appreciable only on a very small scale (see inset). This fine tuning analysis 



9 



shows that each strength grows slowly with 1/e, apart from the transition 1/2^ — > 3/2+, 
whose strength decreases at an even slower rate. In all cases, their smooth behavior allows 
for an extrapolation to asymptotic values through a formula having the same structure as the 
scaling law adopted for the energies (Eq. 5.5). Table II shows that the strengths computed 
at a relatively large e differ very little from the extrapolated ones, which in turn practically 
coincide with the exact values. This rapid convergence is quite significant in view of the 
extreme sensitivity of the transition strengths to even very small components of the wave 
function. 

VI. CONCLUDING REMARKS 

We have shown that the importance sampling is inherent in the iterative algorithm for 
diagonalizing large matrices. The truncation of the dimensions of the matrices it promotes 
is not only kept under strict control, but is also quite severe when the eigenvectors are highly 
localized. For shell model nuclear Hamiltonians, we achieved this localization by adopting a 
correlated basis obtained by partitioning the shell model space into two or more subspaces. 
As clearly illustrated by the calculations carried out on some typical nuclei, the sampling 
so implemented allows to reduce the sizes of the Hamiltonian matrix by at least an order 
of magnitude with no detriment of the accuracy. We get in fact very accurate energies, 
wave functions and E2 reduced strengths. Moreover, it generates extrapolation laws to 
asymptotic eigenvalues and E2 transition probabilities which coincide practically with the 
exact corresponding quantities, whenever these are available. 

It is important to point out that the method is specially effective when applied to 133 Xe, 
having a neutron excess. On the ground of this result, we feel confident that the sampling, 
combined with the use of the correlated basis, will enable us to face successfully the eigen- 
value problem in heavier nuclei, all having a neutron excess. We also like to stress that the 
partition method is specially suitable for enlarging the standard shell model valence space 
through the inclusion of n-particle n-hole correlated configurations. We plan to make such 
an implementation in order to study the intruder states in light as well as heavy nuclei. 
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FIGURES 



FIG. 1. Importance sampling parameters versus the dimensions n of the truncated matrices 
resulting from the sampling. 

FIG. 2. Energy distribution of the J n = + basis states in 48 Cr obtained by the multipartition 
method (MP) and corresponding running sum. This is compared with the running sum of the 
standard shell model (SM) basis states. 

FIG. 3. low- lying J 71 " = 2 + eigenvalues of 108 Sn versus the dimensions n of the truncated 
matrices resulting from the sampling. 

FIG. 4. The same as Fig. 3 for 48 Cr 



FIG. 5. The same as Fig. 3 for 133 Xe. 

FIG. 6. Eigenvalues obtained by the sampling and by energy truncation versus the truncated 
dimensions n. 

FIG. 7. Overlap of the lowest five J 71 " = 2 + and J w = + sampled eigenfunctions with the 
corresponding exact ones in 108 Sn and 48 Cr, respectively. 

FIG. 8. Sampling of the strengths of J w = 2+ -> J 71 = 0+ E2 transitions in 108 Sn and 48 Cr. 

FIG. 9. Sampling of the strengths of some El transitions in 133 Xe. 
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TABLES 



TABLE I. Approximate, extrapolated and exact (when available) energies (in MeV).The ex- 
trapolation parameters are also given. 



A X Jf N Eq b c n E^ Ef xtr) Ej(ex) 

108 Sn 2+ 17467 -3.145 0.0028 0.0032 1867 -3.120 -3.143(001) -3.151 

2+ -2.466 0.0029 0.0034 -2.439 -2.463(001) -2.469 

2+ -2.244 0.0030 0.0021 -2.216 -2.244(001) -2.266 

48 Cr Of 41355 -32.972 0.01721 0.00456 5739 -32.851 -32.955(019) -32.954 

0+ -28.638 0.03450 0.01486 -28.414 -28.604(011) -28.565 

0+ -27.166 0.06038 0.01280 -26.819 -27.166(042) -27.158 

2+ 182421 -32.164 0.019 0.008 14642 -31.952 -32.145(020) -32.148 a 

2+ -29.133 0.0114 0.0009 -29.012 -29.122(016) 

2+_ -28.745 0.0238 0.0037 -28.552 -28.721(053) 

133 Xe l727 125756 -11.280 0.0014 0.0018 6763 -11.256 -11.279(002) 

1/2+ -10.319 0.0019 0.0019 -10.285 -10.317(001) 

1/2+ -9.498 0.0021 0.0019 -9.464 -9.496(007) 

3/2+ 242308 -11.491 0.0009 0.0012 7417 -11.466 -11.490(002) 

3/2+ -10.491 0.0012 0.0012 -10.455 -10.490(003) 

3/2+ -10.292 0.0013 0.0013 -10.255 -10.292(002) 

11/27 497666 -11.222 0.0004 0.0005 6472 -11.192 -11.222(003) 

11/2^ -10.013 0.0006 0.0005 -9.971 -10.012 (003) 

11/27 -9.439 0.0008 0.0002 -9.386 -9.438 (008) 



a from Ref. [21]. 



TABLE II. Approximate, extrapolated and exact (when available) B(E2) (in e 2 fm 4 ). 



A x 


J? 


- J 7 


e 


m 


n / 


B{E2\e) 


B(E2)( extr ) 


B(E2)( ex - ) 


108 Sn 




-o + 


l.xl0~ 4 


1867 


1034 


42.45 


42.68(02) 


42.71 


48 Cr 


2 + 


->0+ 


4.xl0~ 5 


14642 


8144 


226.5 


227.8(8) 


228 


133 Xe 


1/2+ 
V2 + 


-3/2+ 
-3/2+ 


7.x 10~ 5 


6763 


7417 


110.1 
230.9 


109.3(7) 
233.0(3) 
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